Abstract 



In this paper we present a generalization of the Goraychev-Chap- 
lygin integrable case on a bundle of Poisson brackets, and on Sokolov 
terms in his new integrable case of Kirchhoff equations. We also 
present a new analogous integrable case for the quaternion form of 
rigid body dynamics' equations. This form of equations is recently 
developed and we can use it for the description of rigid body motions 
in specific force fields, and for the study of different problems of quan- 
tum mechanics. In addition we present new invariant relations in the 
considered problems. 
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GENERALIZATION OF THE 
GORAYCHEV-CHAPLYGIN CASE 



A. V. BORISOV, I. S. MAMAEV 



In this paper we present a generalization of the Goraychev-Chaplygin 
integrable case on a bundle of Poisson brackets including (co) algebras so(4) 
and so(3, 1), and on the quaternion form of Euler-Poisson equations. Note 
that the generalization on the bundle is connected with the introduction of 
variables on the bundle analogous to Andoyer-Deprit variables. These vari- 
ables are separating for all members of the bundle. We also obtain the L — A- 
pair of the generalized Goraychev-Chaplygin case of quaternion equations. 
In a particular case this pair is reduced to the L — A-pair of the classical 
Goraychev-Chaplygin case. 



Let us consider a generalization of the classical Goraychev-Chaplygin case 
of Euler-Poisson equations 

„> dH dH dH 

M = M x — h 7 x — — , 7 = 7 x 



DM ' 97' ' ' dM 

for the zero value of the area integral (M, 7) = 0. We add a gyrostatic 
moment and singular term to the Hamiltonian 



H = ^{Ml + Ml + 4M 3 2 ) + AM 3 + /i 7 i + i 1 ) 

2 7 3 



where A, /x, a = const. The additional integral is of order three with respect 
to the moments 

F = (m 3 + ^) [Ml + Mj + -|) - iiM ll3 . 

D. N. Goraychev himself in the paper showed the generalization ([|) on 
zero value surface of gyrostatic moment A = (for A 7^ 0, a = it was 
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shown by L. N. Sretensky 0). The complete form of generalization (|T]) was 
considered by I. V. Komarov and V. B. Kuznetsov ||,. They also presented 
some quantum mechanical interpretation of the singular term. 

The Poisson bracket of variables M, 7 is defined by algebra e(3) and has 
the following form 

{M h Mj} = -s ijk M k , {Mi,jj} = -s ijk j k , {7,7,} = 0. (2) 

Let us present the generalization of the Goraychev-Chaplygin case on the 
bundle of brackets of the form 

{M h Mj} = -e ijk M k , {Mi,jj} = -e ijk j k , {7^ 7,} = -xe ijk M k . (3) 

At x — 1 these commutation relations correspond to the algebra so(4). 
We can present the Casimir functions of bracket (|||) in the form 

F 2=1 2 + xM 2 , F 1 = (M, 1 ). (4) 

Now let us construct the variables on bracket ||) analogous to the Andoyer- 
Deprit variables 0. 

Suppose the component M 3 is equal to the momentum 

L = M 3 . (5) 

The variable / canonically conjugate to the variable L ({I, L} = 1) on sub- 
algebra so(3) with the generators Mi, M 2 , M 3 is constructed by integration 
of the Hamiltonian flow with Hamiltonian function H. = L 

dM 1 t , dM 2 

—L = { Ml , L} = M 2 , —1 = {M 2 , L} = -M x , 

Hence using the commutation relation {M 2 , M 2 } = —M 3 we obtain 



M 1 = VG 2 - L 2 sin I, M 2 = VG 2 - L 2 cos I, (7) 

where G 2 = M 2 + M| + Mf is the Casimir function of subalgebra so(3). 

Suppose G is the second momentum and construct the canonically con- 
jugate variable g. We chose H = G as a new Hamiltonian, and the corre- 
sponding flow on the whole bundle C x has the form 

dM c?7 1 
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where g is the variable canonically conjugate to G. 

According to (||), M does not depend on g, and using equations (||) and 
Casimir functions (|) we obtain for 7 the equations 

H a 

7 = + -(Mxe 3 sin 9 + GM x (M x e 3 ) cos 5), 

^(2 H 2 (^) 

a 2 = ° 2 ~ G X 2_L2^ > e 3 = (0. 0. !)> 

where c 2 = xAi" 2 + 7 2 , and if is the traditional notation of the area inte- 
gral H = (M,7) = c x . 

Thus (|), ([?]), (|9|) define the symplectic variables on the whole bundle C x 
corresponding at x = 0, c = 1 to the well-known Andoyer-Deprit variables 
in rigid body dynamics. 

Using (|) , (H) we obtain the generalization of the particular Goraychev- 
Chaplygin integrable case on the bundle C x . We chose the Hamiltonian in 
the form 

H = ]-(G 2 + 3L 2 ) + XL + a(cos/cosg + sin Z sing), (10) 

2 Gr 

where a, A are constants. 

In ( |i~0| ) we add the linear with respect to L term. It is interpreted on 
algebra e(3) as the component of the gyrostatic moment ||. 

We can separate variables in system (|i~0|). Indeed, let us apply the canon- 
ical change of variables 

L=p 1 +p 2 , G = pi-p2, qi = l + 9, 02 = 1- 9- (11) 
Now Hamiltonian ( |l0| ) can be presented in the form 

H = \ P ^^ - + -^( Pl sm gi + P2 sing 2 ). (12) 

2 Pi - P2 Pi - P2 Pi ~ P2 

Using (0), (|^), © we present Hamiltonian ( |10|) as a function of variables M, 7 
on the zero surface of area integral (M, 7) = H = 

H = i(M 2 + M 2 + 4M 3 2 )+AM 3 + /^. (13) 
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The additional integral in this case has the form 

F = (m 3 + {Ml + Ml) - fiM,^- . (14) 

On algebra e(3) we have [— y | = 1 and obtain the classical Goraychev- 
Chaplygin integrable case. In the classical case this method of variables' 
separation was suggested by V. B. Kozlov 0. 

Singular term (|I|) is generalized on the bundle in the following way: 

H = -{Ml + Ml + AMI) + AM 3 + /i^ + - a7 ~ 



r --, ■ < ■ , " ' '-| 7 | ' 2 7 f ' 

F = (m 3 + ^) (m 2 + M 2 2 + a^§) - //Mi 



'173, 



7s- 

although variables' change (|TT1) in this case doesn't produce the separation 
of variables (at least we don't know such separation). 

Recently Sokolov and Tsyganov present a new generalization of particu- 
lar integrable case (||) on bracket (||) with the Hamiltonian containing the 
quadratic cross terms with respect to M , 7. 

The most general form of the integrable family in this case is presented 
as the following Hamiltonian 

H = UmI + Ml + AMI + 4) + AM 3 + A*i7i + A<272 + , . 
2V " 7 3 y (15) 

+ ai {2M 3ll - M173) + a 2 (2M 3 7 2 - M 273 ). 
The additional integral is 

F = (m 3 + ai 7l + «272 + ^) (M 2 + Ml + -|) - (/iiMi + /i 2 M 2 )7 3 . (16) 

At e = we can separate variables in Hamiltonian (fl5|) using the linear 
combination of Andoyer-Deprit variables analogous to (pi]). Indeed, we can 
show that 

H = — - — (/(pi,9i) -#(P2,<? 2 )), 
P1-P2 

f = 2aip\ + 2 + ai sin gi + a 2 cos q^\p\ + (yUi sin g x + /i 2 cos <?i)pi, (17) 
g = 2aipl + 2 ^ - ai sin g 2 - a 2 cos q 2 ^jpj - (/ii sin q 2 + fi 2 cos q 2 )p 2 . 
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As we note above, we don't know such separation at e ^ 0. 

If the gravity field is absent (/xi — f/,2 — 0), then the integral is presented 
product of factors F = k x k 2l where 

h = a M 3 + a lll + a 2l2 + ^, k 2 = + M% + ^. (18) 

^ 73 

We can easily show that the equations of motion for ki, k 2 have the form 
fci = 2(ai7 2 - a 2 ii)k 1 , h = -2(a^ 2 - a 2 ^i)k 2 . 

(M,-y)=0 

Thus the equations ki = and k 2 = define the invariant relations of 
system (|15|) at //i = //2 = (and fc 2 = is the particular invariant relation 
on zero level of (M, 7) = 0, and fci = is general invariant relation). It seems 
that, invariant relations (|I8"D and the additional integral in the form F = k\k 2 
where introduced by authors. 



Using the observation from the first section we can generalize integrable 
case ( pT5| ) and invariant relations on bundle of brackets (|3|). We shall make 
the following substitution in Hamiltonian (|15|), integral (|16|), and invariant 
relations (|i"8|) 

n o, ItI = \Jil + il + il- 

In this case at e = we can also separate variables using the analog of 
Andoyer-Deprit variables and substitution (|Tl|) . 

We should note that unlike the generalization of Kovalevskaya case |l|, 
where the Hamiltonian and integral explicitly depend on the parameter of 
bundle x, in the generalization of Goraychev-Chaplygin case the integrals 
don't depend on the parameter of bundle. 



Let us consider another possibly not so natural case of motion equations of 
rigid body with a potential linear with respect to Rodrig-Hamilton parame- 
ters and not with respect to direction cosines |l| 

1 3 
H = -(AM, M) + ^VjAj, n = const, (19) 

8=0 
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The equations of motion in this case are 



M = M x h -A x — + -- — A - —An -7; — , 

dM 2 9A 2 3A 2 dA ' 

A °-"21 A '9mJ' A " 2 X 9M + 2 Ao 9M' 

where A = (Ai, A 2 , A3). 

These equations are Hamiltonian one with the Poisson bracket 

{M h Mj} = -e ijk M k , {Mi, A } = -A i; 
{Mi,\i} = --(eijkXk + SijXo), i,j,k = 1,2,3, 



(20) 



(21) 



corresponding to the Lie algebra so(3) ® s R 4 . The relations between quater- 
nions and Euler angles have the form 

9 ijj + if .9 - f 

An = cos - cos , Ai = sin - cos , 

u 2 2 2 2 

. 6 . if; — if . 9 . ip + f 

Xo = sin - sin , A3 = cos - sin . 

2 2 2 2 

In the classical mechanics such potential aren't presented because its depen- 
dance on the body position is nonunique (more exactly it is double-valued 
function). To justify the study of such equations we can refer to the problems 
of quantum mechanics and point masses dynamics in curved space S 3 , and on 
some formal technics of L — A-pair construction [||, ||. It turns out that after 
the reduction of order of system ( |l9| ) we obtain the general Euler-Poisson 
equations with additional terms having different physical interpretations QTJJ . 

As an interesting feature of system (|19|) we note that using linear with 
respect to A« transformations we can reduce the general form of the potential 



V = J2 r ^i ( 22 ) 



to the form 



V = r A . (23) 
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Indeed, the linear transformations of quaternion space Aj (preserving com- 
mutation relations and quaternion norm) of the form 



A = R~\r \ + nAi + r 2 A 2 + r 3 A 3 ), 

Ai = i2 -1 (r Q Ai - n\ - r 2 A 3 + r 3 A 2 ), 

A 2 = R~ 1 (r Q X 2 + r x A 3 - r 2 A - r 3 Ai), (24) 

A 3 = i2 -1 (r Q A3 - ^iA 2 + r 2 Ai - r 3 A ), 
R 2 =rl + r\ + r 2 2 + r 2 3 

reduce potential (p2|) to form (p3|). The existence of such linear transforma- 
tion is the remarkable property of quaternion variables and bracket (|2ll) , the 



analogous transformation does not exists for the brackets of algebras e(3) 
and so(4). 

In general dynamically unsymmetrical case a\ ^ a 2 7^ a 3 7^ ai system (|T9|) 
seems to be nonintegrable and none of two necessary additional integrals 
exists. However, this was not proved, and the proof by various reasons is 
not simple. Note that even the application of the Kovalevskaya method for 
system (|i"9|) is not quite analogous to the classical Euler-Poisson problem. 
Generally speaking even in the Euler-Poinsot case the solution branches on 
the complex plane of time (with the exponent 1/2). 

At ai = a 2 the linear integral always exists 



F x = M 3 (r 2 + r\ + r\ + r\) + N z {r\ + r\ - r 2 - r|) + 
+ 2N 2 (r 1 r - r 3 r 2 ) - 2N 1 (r 1 r 2 - r r 3 ), 



(25) 



where Ni are projections of kinetic moment on the fixed axes. For poten- 
tial ( |23|) this integral has the natural form 

Ft = M 3 - N 3 . (26) 



It turns out that (linear) integral (|26|) corresponds to the cyclic variable (f + 
ip Rauss reduction with respect to this cyclic variable results in Hamil- 
tonian system on algebra e(3) on zero surface of area integral (M",7) = 
with the Hamiltonian 

H = \{Ml + Ml + a 3 M 2 ) + c(a 3 - 1)M 3 + r o7 2 + (27) 
z 1 7 3 
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where c is a constant value of integral (|26|). Hamiltonian (|27[) corresponds to 



the addition of a gyrostatic term linear with respect to M and singular term 
into the general Euler-Poisson equations. The physical meaning of singular 
term is discussed in [l]]. 

Integrable cases of system (|i"9|) (equivalent to integrable cases of sys- 



tem (|27D ) are presented in 0. Here we just present the generalization of 
GoraychevChaplygin case. 

Hamiltonian and additional integral are 

H = l -{Ml + M 2 2 + AMI) + n)A , 

F 2 = M 3 (M 1 2 + Ml) + r (M 2 Ai - MiA 2 ). 

Integral F 2 commute with integral Fi. Under the reduction to system fl27j) 



this case becomes a member of generalized family ([[]). 

Remark 1. Adding a constant gyrostatic moment along the axis of dynam- 
ical symmetry in (p8| ) we obtain an integrable case corresponding to the 
generalized Sretensky case in Euler-Poisson equations. Additional integrals 
can be easily obtained with the help of the lifting procedure described in ^j. 

Note that the indicated "Goraychev-Chaplygin case" for the quaternion 
Euler-Poisson equations is the general integrable easel Thus we can use 
it for some algebraic constructions (L — A-pair constructions and all that) 
described below. 



Lax representation of systems on algebra su(2, 1). We start from 
some formal method of L-A pair construction based on existence of two 
compatible Poisson brackets for the same system of Hamiltonian equations 
on Lie algebra; in this case the system is called bihamiltonian. The detailed 
exposition of this method and its applications to various problems of rigid 
body dynamics are presented in P, |J. 
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Let's consider space C of complex matrixes 3x3 with the basis 



Mi = 


(0 \i 











M, 



Pi 



P 3 = 



(-¥ o 




















\-5* 















\±a;i 





M 2 = 













IVL 



P 4 = 










V 













f 










^-5^ 









\0 \x% 





(29) 



With respect to the standard matrix commutator [• •] such matrices gen- 
erate a semi-simple algebra with Cartan decomposition C = H + V, where 
subalgebra H = su(2) © su(l) is generated by matrices Mj, and V = C 2 is 
generated by matrices Pj. Here a; is a parameter determining some bundle 
of algebras linearly dependent on x. For x > these algebras are isomorphic 
to algebra su(3), for x < they are isomorphic to algebra su(2, 1), at x — 
they are isomorphic to the semidirect sum (so(2) © su(l)) © s C 2 . 

Using the semi-simplicity we can identify algebra with coalgebra by means 
of the inner product (Killing form) 



4 = -Tr(X-Y) 



X, Y 6 £. 



(30) 



Let's denote coordinates in coalgebra as mi, m 2 , m 3 , m 4 , pi, p 2 , p 4 , 
then after the identification we obtain a matrix (an element of algebra): 



X 



'-i(m 3 + m 4 ) im 1 + m 2 ~ 
im 1 -m 2 i(m 3 -m 4 ) -{p 3 -ipi) 
-pi - ip 2 -p 3 - ip 4 2im 4 



The corresponding Lie-Poisson bracket, more precisely, the bundle of 
brackets linearly dependent on parameter x has the following form for the 
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coordinate functions of coalgebra 



{mi, rrij} = e ijk m k , {m u m 4 } = 0, 

1 1 

{m h pj} = -(e ijk p k - SijPt), {m u p 4 } = -pi, 

{Pi, Pj} = \x{e ijk m k + 3e ij3 m 4 ), {p h p 4 } = -\z 




In correspondence with the general method developed in || |j, we shall 
construct L matrix, and its invariants define a commutative set for the whole 
family of brackets { }<?+A({ }a + { } a ), where a e V is a shift of the argument. 
The bracket {}e is different from ([H]); in this bracket variables pi are pairwise 
commuting. Let's assume x — 1 restricting ourselves to the problem with a 
real dynamical sense. Then we can present L-matrix as: 



Among invariants of matrix L under arbitrary shift there is a linear on rrii 
function of the form 



L = (hA + v + aA 2 ) 



where 




(32) 



Fx = m 4 a 2 + (m, j a ), 



(33) 



4 

where a 2 = 2, fl2 ; m — ( m i ? m 2, m s), and the components of vector *y a are 



i=l 



the functions of of the form 



7 a = (2(aia 3 + a 2 a 4 ), 2(— a 2 a 3 + a\a 4 ), a 3 + a 4 — a x — a 2 ). 
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Let's suppose ai = a 2 = 0. In this case integral (|33|) has the form 

F\ = m 3 + m 4 . (34) 

The following square-law invariant of matrix L, we will choose as a Hamil- 
tonian 

F 2 = Tr(h 2 + 2va) = m\ + m\ + m\ + 3m 2 + a 4 p 4 + a 3 p 3 . (35) 

It defines some (formal) integrable Hamiltonian system on the family of 
brackets {}g + A({}a + {} a )- We can easily find matrix A for this system ||. 



The reduced system and nonlinear Poisson structure. In order to 
proceed from the found formal system to the generalization of Goryachev- 
Chaplygin case |7|], we shall make the reduction using linear integral (Q) 
directly in obtained L and A matrices. In the general case this reduction 
is not Poisson reduction ||. We shall make the following substitution in 
matrix L fl30|) and in Hamiltonian (|35"D 



m 4 = — m 3 + c, c = const . 

We obtain L matrix and Hamiltonian of integrable system on subalgebra m l5 
m 2, rm, Pi, P2, P3, P4 with gyrostat, the gyrostatic moment being equal to c: 

— iXc (imi+m 2 )A Pi—ipi 

(im 1 -m 2 )A i\(2m 3 — c) p 3 — zp 4 +(a 3 +ia 4 )A 2 | (36) 

-Pi-ip 2 -p 3 -ip 4 -(a 3 +a 4 )A 2 -2i(m 3 -c)A 



H = m l + m 2 + 4m 3 — 6m 3 c + 2a 4 p 4 + 2a 3 p 3 (37) 

(i(4m3 — 3c) — %ra\ — 1112 
—im\ + ni2 —i{Am-i — 3c) 
as + ia^ 

Let's consider Poisson structure, defined by a bracket {}g. The corresponding 
Hamiltonian system with Hamiltonian fl37|) has another linear integral 




F 3 = m 3 - (m, 7) (38) 

7 = (2(p 2 P4 + pm), 2(p 2 ps - pm), pI+pI-pI- pD), 
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using this integral we can carry out a usual procedure of order reduction (such 
as Routh reduction or reduction on moment). We can carry out the reduction 
in the simpler way using algebraic form, if we choose as new variables the 
integrals of vector field v = {•, F 3 } of the form || 

M lPl + M 2 p 2 M 2Pl - M lP2 

K i = / 5 5 , K 2 = K 3 = M 3 , 

V Pi + Pi VPi + Pi (39) 



having the nonlinear commutation. We describe the transformation 



in 



-^4 

b 3 (,4UJ 



where F 3 = (X, s)s 3 is a Casimir function of structure (|40|) . 

At zero value of "the area integral", i. e. for F 3 = (K, s)s 3 = bracket ( fHf ) 
coincides with a usual LiePoisson bracket, defined by algebra e(3) = so(3) © s M 3 , 



and Hamiltonian (|37|) in variables fl39|) coincides with the Hamiltonian of the 
classical Goryachev-Chaplygin case [0] 

#* = ~# = i(i^f + iT| +4Kf) - a 4 s 2 - a 3 si, a 3 , a 4 = const . (41) 

Goryachev case with a singular term. We can remove the nonlinearity 
of structure fllCf ), appearing for F 3 ^ on the fixed level F 3 = c using the 
transformation 

L = K — c— . 

S3 



After the transformation bracket (|40| ) gets the form of algebra e(3), and 
Hamiltonian if* gets the form 

H* = \ (Lj + L\ + ALl) ~ «4S 2 - a 3Sl + 3cL 3 + ~ . (42) 

Hamiltonian ([42] ) can be interpreted as some generalization of the Goryachev- 
Chaplygin case, at (L, s) = 0, with the additional terms, linear on L 3 and 
corresponding to the gyrostatic moment. The integrable generalization with 
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the gyrostatic moment only was described by L. N. Sretensky in ||, the gen- 
eralization with singular potential only was presented by D. N. Goryachev 
himself |7j, the general case, when it is possible to add to the Hamilto- 
nian both terms with arbitrary independent coefficients, was described in 
paper |g. 

Thus, the presented L-A pair is valid for the generalizations of the 
Goryachev-Chaplygin case. It is different from the mysterious L-A pair 
described in paper || which is obtained by deletion of a row and column 
from the relevant pair of the Kovalevskaya case. 

The authors thank A. V. Bolsinov, A. V. Tsyganov, V. V. Sokolov for the 
valuable remarks and numerous discussions. 
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